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1 The iterative calibration routine uses trial 
and error to compare the set of 4 measured 
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2 In the dual -wedge scanner coordinate system, 8 

the input vector S is coll inear with the 

x axis, which is the axis of rotation for the 
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3 Using simple approximation, counterrotating 9 

two equal wedges of the scanner produces a 

straight line on the z axis. 

4 A closer look at counterrotating equal wedges 12 

shows errors of 30 meters at 10,000 meters 
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greater distortions to the scan generated. 

6 Rotating the second wedge independently pro- 16 

duces the ellipse shown. The desired point 

can then be arrived at by rotating the wedges 
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I. INTRODUCTION 

The ability to precisely position a laser beam along one of 
many desired lines of sight is vital to many laser radar systems. 
It is often beneficial to house the equipment in a small package 
as well, so that, for example, the dimensions of an output window 
which the scanned beam passes through may be made as small as 
possible. When limited spatial coverage is sufficient and this 
packaging constraint exists, a scanner using two wedges which may 
be rotated independently about the direction of the input beam is 
an attractive choice. A 30.5-cm-diameter germanium dual-wedge 
scanner has been built for a CC^ doppler lidar used on NASA's 
CV-990 aircraft for meteorological research 1 and was tested on 
that aircraft during a flight test in 1981. 2 The success of this 
test demonstrated that precise, reproducible pointing could be 
achieved with this device. Two problems were observed during this 
test and the subsequent analysis: 

(a) The dependence of the deviation on the angle of incidence 
at the second wedge, while small, was sufficient to produce an 
error greater than that caused by other sources, and more accurate 
positioning calculations were therefore required. 

(b) A more accurate alignment procedure was required to 
establish the relationship between the scanner axes and a known 
reference frame. 

Since the desired reference frame was determined by the inertial 
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navigation system (INS) of the aircraft and no reference surfaces 
were available near the scanner, it became obvious that this 
alignment must be performed relative to the INS-indicated air- 
craft orientation. 

If the scanner alignment parameters— wedge angles, wedge ori- 
entations, and the axes of the scanner and input beam direction- 
are known exactly, the line of sight of an output beam can be 
determined from the indicated wedge positions as shown in 
Figure 1(a). This calculation makes use of Snell's law in vector 
form 3 and coordinate rotation matrices. More frequently, it is 
desired to position the wedges to produce a specified line of 
sight as shown in Figure 1(b). This can be accomplished with 
only slightly greater difficulty as will be shown later. 

Before either of these calculations can be performed, it is 
necessary to determine the alignment of the scanner. In prin- 
ciple, if enough combinations of measured lines of sight and 
indicated wedge positions are determined, a set of simultaneous 
equations can be generated and solved for all the unknowns among 
the alignment data as shown in Figure 1(c). This is accomplished 
by locating the beam on retroreflector targets at surveyed points 
and recording the indicated wedge positions which produced each 
known line of sight. In practice, the inversion process is 
difficult, and is solved by an iterative approach as shown in 
Figure 1(d). 
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(a) Calculating the Line of Sight 

Figure 1. The iterative calibration routine uses trial and 
error to compare the set of 4 measured data points 
with trial data points calculated with aircraft 
orientation, scanner setting, and scanner alignment 
data. 


3 










(b) Setting the Scanner to a Line nf Sight 

Figure 1. The iterative calibration routine uses trial and 
error to compare the set of 4 measured data points 
with trial data points calculated with aircraft 
orientation, scanner setting, and scanner alignment 
data. 
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(c) Calibrating the Scanner 

Figure 1. The iterative calibration routine uses trial and 
error to compare the set of 4 measured data points 
with trial data points calculated with aircraft 
orientation, scanner setting, and scanner alignment 
data. 
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(d) Iterative Calibration Routine 

Figure 1. The iterative calibration routine uses trial and 
error to compare the set of 4 measured data points 
with trial data points calculated with aircraft 
orientation, scanner setting, and scanner alignment 
data. 
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The approach to precision pointing will be described beginning 
with the generation of a simple line scan, followed by a descrip- 
tion of the full-coverage pointing algorithm. Next, the iterative 
calibration algorithm will be described along with its applica- 
tion to a model scanner. Finally, the calibration algorithm will 
be applied to a limited set of calibration data collected during 
preparation for the 1981 flight tests, and the results compared 
to actual flight test results. 

2. LINE SCANS 

One simple application of a dual -wedge scanner is the genera- 
tion of straight line scans. As shown in Figure 2, this scanner 
uses two wedges which have wedge angles that are as nearly equal 
as possible arranged with their flat sides together. The first 
wedge produces a certain deflection of the beam in the plane of 
incidence. As the wedge rotates about the beam direction, the 
plane of incidence, and thus the plane in which the output beam 
emerges, rotates as well. If the x axis is defined as coll inear 
to the beam and is thus the axis of rotation as the first wedge 
is rotated through 360°, a point at a specified range along the 
deflected beam traverses a circle parallel to the y-z plane 
centered on the x axis (Figure 3, box 1). With this wedge at 90°, 
the second wedge, W 2 , is introduced. If the deviation were 
independent on the angle of incidence, a rotation of W 2 through 
360° would produce a second circle centered on the first circle 
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LASER 


Figure 2 



In the dual -wedge scanner coordinate system, 
the input vector S is coll inear with the 
x axis, which is the axis of rotation for 
the two wedges. 


8 




) 

r 

f 

SV 

S .. B 

■ns 


o 

o 

II 

CD 

<J 

vs 

0 

A# = 90° 

0 

■ 



■ 

J V 

oo 

II 

CD 

<] 

0 

II 

CD 

<] 

0 


/\ Q 2 0 l 

Hi 

A r 2 0 l 

vs 

A0= 360° 

J 

© 

9 1 ' 

V 

Ad - 720° 

7 )fV 
© 


Figure 3. Using simple approximation, counterrotating 
two equal wedges of the scanner produces a 
straight line on the z axis. 
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at 0j = 90°. For moderate wedge angles, this is approximately 

true, and was, in fact, assumed during the 1981 tests. 

With 0^ = 90° and 0 ^ = 90°, the first point of the scan would 

be at maximum z or twice the radius of the circles out from the 

origin. Next, the wedges would be counterrotated' 45° each, 0^ 

counterclockwise and 0^ clockwise looking along the direction 

of the impact beam until A0 = 90° (box 2), and the beam would be 

at point 2 on the z axis. At A0 = 180°, point 3, which coincides 

with the origin, is reached, and further counterrotation moves 

the beam to the minimum z at point 5. Another 360° results in 

complete counterrotation and a return to the maximum z value. In 

this way, a straight line scan of length equal to twice the 

diameter of the circles would be produced. 


However, while the angle of incidence at is constant for 
all rotation angles, it is not at since the beam is refracted 
by W^ in one of its many possible orientations. For this reason, 
the beam path in this scan is not precisely straight. The exact 
form can be obtained using the vector form of Snell's law, 3 


Vi ■ r jV + |VH [hvv*] - vVM "j • (i > 

where at surface j, the vector S, is a unit vector parallel to the 
incident beam, S. + j is parallel to the refracted beam, N is nor- 
mal to the surface, and r^. = N ^ /N ^ is the ratio of the index 

of refraction on the side with the incident beam to that on the 


side with the refracted beam. If the inner wedge surfaces are 
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parallel, they may be neglected; therefore, the two surfaces to 
be considered are the inner surfaces of the inner wedge and the 
outer surface of the outer wedge, described by 


A A A /V 

N. = cos w. i + sin w. (cos 8. j + sin 0. k), (2) 

J J J J J 

with a beam incident along the x axis described by 

A A 

s x = 1. (3) 

For equal wedges, the scan pattern follows the bow-tie curve 
shown in Figure 4. This curve, drawn with w^^ = w 2 = 3.3275°, 
shows that the errors in the y direction can equal approximately 
30 meters at a range of 10,000 meters. This error reduces to 
zero whenever 9^ and 0 2 are each at 90° intervals. Figure 5 
shows the effects of unequal wedges on the curves and illustrates 
some of the problems associated with precise positioning using 
a dual -wedge scanner. 

3. FULL COVERAGE PRECISE POINTING 

The problem of precise pointing could be solved by inverting 
the vector equations for refraction at two surfaces to solve for 
0j and 0 2 , but this is a difficult or impossible task. A realis- 
tic approach has been developed using a two-step process which 
may be described in terms of its deviation angle from the x axis 
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Figure 4. A closer look at counterrotating equal wedges 
shows errors of 30 meters at 10,000 meters 
for a scan 8,000 meters wide. 
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Figure 5. Counterrotating unequal wedges produces 

greater distortions to the scan generated. 
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and its rotation about that axis. The line scan shown earlier 
may be used to generate any deviation given by 

dev = cos 

The key to the algorithm is that the equation for the deviation 
as a function of ~ 9 i can inverted. If one of the wedge 

angles is known, this inversion also provides a rotation angle. 
The second step is to rotate both wedges by the additional amount 
required to achieve the desired rotation. 

The two steps are defined precisely as follows. It is con- 

A 

venient to assume that 0, = 0° initially. Then the beam S, = 

As As as i 1 

li + Oj + Ok is refracted through Wj, whose unit normal is given 

A A A 

by N, = (cos w.) i + (sin w, ) j. The vector equation for a 

111 As As 

refracted beam is given by equation (1). Substituting S^.N^ = 
cos w 1 , the beam refracted through the first wedge is described 
by: 



A 



r.i +Vl-nf (sin^ w.) - r. cos w. 


£(cos w^) i + (sin w^) jj. 


(4) 


As expected, there is no z component, since 0^ =0°. 

A 

To derive an equation for the output vector, and solve 
that equation for the rotation angles of the wedges, equation (4) 
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is used again in equation (1), along with N 2 , calculated using 
$2 = A0, the difference angle. This, in effect, is the same as 
counterrotating the two wedges except that the frame of reference 

A 

rotates with 9^ (see Figure 6). N 2 is given as 

A A A A 

N 2 = (cos w 2 ) i + (sin w 2 cos A0) j + (sin w 2 sin A0) k. (5) 

Since, by definition, the final unit vector along the beam from 
the second wedge will have an x component x^ = cos (dev). 


cos (dev) = r 2 x 2 


( 6 ) 


+ 


{yi-rf 


r 2 £ 2 i 2 j 


COS W 2 , 


which is related to the angle A0 only by its dependence on S 2 .N 2> 
This dot product is obtained using the two-step process shown 
below. Defining the quantity in brackets in equation (6) as 


cos (dev) - r,x 


cos w 2 


C p 

1 — T r~~ ' 

2 

cos (dev) - r 2 r^sin (w^) + cos (w^)* 1 

^1-r 2 sin 2 w 1 

cos w 2 


( 7 ) 


A /\ 

and solving equation (6) for S 2 .N 2 , 
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e 1 = 90° 


9, = 90° -A0 
c 0° < A0 < 360° 


.(3) DESIRED 
POINT 



Y AXIS 


Figure 6. Rotating the second wedge independently 
produces the ellipse shown. The desired 
point can then be arrived at by rotating 
the wedges simultaneously. 
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S 2 .N 2 


1-r 2 - 6 1 

2<$r 0 


Recalling equation (5) for N 2 , the difference angle is given by 

Jl 


“a a 


A8 = COS 


_ 1 j [_S 2 .N 2 - x 2 cos w 2 
y 9 sin w 9 


( 8 ) 


Since the initial condition was 0^ = 0°, the beam has been 
rotated by an amount equal to the angle between the final z and 
y components, 9-j. = tan”* (z^/y^); the desired rotation, 8, may be 
achieved by setting = 0 - 0j and 0 2 = 9 l + A0 ’ 

4. ALIGNMENT AND CALIBRATION 

The initial application of this system was to position a dop- 
pler lidar beam in order to measure wind fields from an aircraft. 1 * 
On an aircraft, there are several reference frames which may be 
useful. Transforming from one frame to another involves only the 
rotation about the three primary axes of the aircraft itself. 

These rotations, shown in Figure 7, are defined individually. 

The heading is a negative rotation (using the right-hand rule) 
about the y axis; heading is defined by the convenient compass 
angles referenced to true north. The pitch is a negative rota- 
tion about the x axis; a positive pitch indicates an upward move- 
ment of the nose of the plane. The roll is a positive rotation 
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Figure 7. Calculating aircraft alignment involves 
three rotations about the coordinate 
systems of heading, pitch and roll. 
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about the z axis; a positive roll indicates that the right wing 
of the aircraft has been moved down. Scanner alignment is done 
in reference to the aircraft, while points are surveyed in the 
earth reference system. In order to perform the necessary cal- 
culation for alignment, it is necessary to convert the measure- 
ment points from earth to aircraft coordinates. Then, just pre- 
vious to comparison, they must be converted back into earth coor- 
dinates. The earth-to-aircraft process involves two steps: 
first, they must be "unheaded"— rotated positively about the y 
axis, an amount equal to 90° less the measured true heading angle; 
second, the vectors must be "unrolled"— rotated negatively about 
the z axis, an amount equal to the measured roll angle. It should 
be clear to see that reversing these calculations affects the 
vectors to return to earth coordinates. The apparent omission of 
pitch correction is deliberate; the rotation of the wedges 
accounts for this movement. 

There are at least six critical parameters for the correct 
alignment of the scanner: w^, W£, 9^ , 0£ , ROLL and THG (true 
heading). In order to precisely position an output beam, it is 
necessary to align the axis of the scanner and the 0° reference 
marks of both wedges with the three axes of the inertial naviga- 
tion system, as well as to ascertain the wedge angles of both 
wedges. Since it would be difficult to locate the required 
reference surfaces and perform the installation accurately, a 
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method has been developed to determine these angles from measure- 
ments made after the scanner is installed in the aircraft; speci- 
fically, a program has been written which calculates the wedge 
and rotation angles as well as the errors in roll and heading 
using four selected scan vectors at which the beam is positioned. 

It would be possible to use the four measured parameters, 
e. n( j, ROLL^j, PITCH..^, and THG..^ (where the subscript "ind" 

A A 

denotes "indicated"), with the measured vectors through to 
calculate errors between actual and predicted lines of sight, and 
through iteration directly adjust the six parameters. This pro- 
cess has two primary faults: the mathematics are complex, and in 

general, all the vectors are affected by a change in any of the 
six parameters. To alleviate these problems, a set of six linear 
combinations of the parameters was chosen, through Pg, and with 
them six check parameters, CK^ through CKg, so that each check is 
most strongly dependent on its own scan parameter. This process, 
which was accomplished through intuition and trial and error, is 
similar to the diagonal ization of a matrix representing a system 
of linear equations. This corresponding dependence makes it pos- 
sible to adjust the scan parameters, calculate new vectors 
through Sg and then their trial check parameters, and finally com- 
pare these check parameters to those of the actual vectors. In 
this way, each scan parameter is adjusted to reduce its own check 
parameters with a minimum of interaction among the parameters. 
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The two sets are listed below. 


Scan Parameter 
P(l) = |w 2 | + | Wj | 
Total wedge angle 


Check Parameter 

Ck ( 1 ) = cos'^S^ .S^) [signed] 

Distance between two points near 
maximum deviation 


P(2) = |w 2 1 - |w x | 

Difference between 
the two wedge angles 

Ck (2 ) = cos _1 (S 3 .S 4 ) [signed] 

Distance between two points near 
minimum deviation 

P(3) =02-0! 
Rotation difference 

-1 

Ck (3) = — where D! = cos 

A A 1 A A 

(Sj.Sgl-cos -1 (S^ - S 2 ) and 

-1 ^ ^ “1 

D 2 = cos (Sj.S^l-cos 

(s 4 .s 2 > 

Difference in deviation between 
two extremes of the scan 

P (4 ) = + 

Ck(4) = tan'^y^-y^/DIST where 

Rotation sum error 

DIST =y(x 1 -x 2 ) 2 + (z 1 -z 2 ) : 


Height difference between hori- 


zontal extremes of the scan 

P(5) = Actual roll — 

Ck (5) = [cos -1 (yo)+cos _1 (y. )]/2 

Indicated roll 


actual heading 

Mean elevation angle near center 

Roll error 

of scan 

P (6 ) = Actual heading— 

■ Ck (6) = [cos -1 (Xt)+cos -1 (x,)]/2 

Indicated heading 

Heading error 

Mean heading near center of scan 
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These six check parameters are calculated from the four surveyed 
lines of sight, and those values are compared to the check 
parameters of scan vectors calculated from the scan parameter 
values. The errors between the two sets of checks are reduced 
by correcting the scan parameters, one at a time, each using its 
own check, until by successive iteration a final set of scan 
parameters is reached. This process has been simulated at 
Raytheon Company's Scientific Computer Center, and the results 
are shown in Figure 8. 

It was assumed above that the incident beam was parallel to 
the axis of rotation. The case for which this is not true has 
also been investigated, and it was discovered that most of the 
errors created by an initial beam rotation around the y and 
z axes are absorbed later in the scan parameter values of head- 
ing and roll errors, respectively. In other words, the deviation 
produced by the wedges is almost independent of the angle of 
incidence for the small alignment errors investigated. The 
following values were calculated to support these claims. 

No Rotation Rotation 1° Rotation 1° 

Correct P's about y axis about z axis 


Wedge Sum 

6.725 

6.730 

6.725 

Wedge Difference 

-0.04 

-0.035 

-0.04 

Rotation Difference 

1.6 

1.8 

1.6 

Rotation Sum 

0.6 

0.7 

0.6 

Roll Error 

0.33 

0.335 

1.37 

Heading Error 

0.75 

1.7? 

0.75 
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WEDGE SCAN PARAMETERS IN DEGREES 




ITERATIONS 


Figure 8. As the iterations progress, each of the 
parameters approaches its correct value. 
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As shown, the rotation about the y axis causes a corresponding 
increase in the heading error, with slight errors also occurring 
in the other parameters. These additional errors are expected, 
since the points selected for the calibration are all close to 
the x-z plane. The roll error absorbs completely the error from 
the rotation about the z axis. 

Using these two misaligned parameter sets, scatter plots of 
errors (see Figures 9 and 10) showed that position errors at 
10,000 meters of range were less than 16 meters. The vertical 
axis on the plots corresponds to the z axis. The circles show 
the outer limits of the 20° scan; therefore, the points inside 
the circles represent all points reachable by the scanner. Each 
point in a circle is represented by a box containing 0 to 16 lines. 
The number of lines is proportional to the error in positioning 
the scanner at that point. A black box, which contains all 16 
lines, represents the maximum error in meters at 10,000 meters 
for each plot, and is shown in the legend. 

5. 1981 TEST RESULTS 

In 1981, a dual -wedge lidar system was developed for NASA's 
Severe Storms Doppler Lidar Flight Program. This system used 
two germanium wedges to point a C0 2 laser beam out the left side 
of a CV-990 aircraft in order to collect wind field data. This 
was done by collecting 1 ine-of-sight doppler velocity data with 
two sequential sets of pulse trains along lines of sight 20° 
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Figure 9. Results better than 14 and 10 meters at 
10,000 meters are acquired with the two 
misaligned parameter sets if the rotation 
angles are unrounded. 
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Figure 10. Due to the one-tenth degree quantization 
in positioning the scanner, the actual 
errors are 15 and 12 meters at 10,000 
meters. 
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forward and 20° aft of the nominal left-looking line every 
1.2 seconds. Then a grid of wind velocity vectors was created 
by combining the doppler measurements along the fore and aft 
lines at all their intersecting points (see Figure 11). 

Previous to the system's installation, the two wedge angles 
were measured to be = 3.3264° and = 3.3206°. Upon instal- 
lation, the system was tested for alignment by directing the 
beam to surveyed points and collecting the associated parameter 
values (see Figure 12). From this set of points, four were 
chosen to be processed through the alignment program to determine 
the correct parameters. Because the proximity of the two center- 
most points prevented the program from successfully finding the 
wedge-difference parameter P^, this parameter was set to zero 
throughout the run, which is equivalent to assuming equal wedge 
angles. The alignment program arrived at values of = 

3.3275°, which are 0.21% and 0.03% from the values measured in 
the laboratory. 

As noted before, the two innermost vectors were too close to 
each other for best possible results. For the greatest accuracy, 
the two outermost points should be separated by no less than 4.6° 
and the two innermost points by at least 0.35°. Even these non- 
ideal choices have been demonstrated to produce successful 
results in the model tests. 
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Figure 11. Fore and aft looking scans create a grid of 
wind velocity returns in the Severe Storms 
measurement system. 

a V 
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OUTER LIMIT 
OF COVERAGE 



COMMUNICATION 



Figure 12. The acquisition and selection of data points 
is restricted by scanner coverage, physical 
obstructions, and human limitations. 
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6. CONCLUSION 

The Severe Storms tests, which ran from 12 June 1981 to 
30 July 1981, used the dual -wedge system successfully to position 
a laser beam for vector wind velocity measurements. An error in 
the data that was originally attributed to misalignment seems to 
have been caused at least in part by the straight line deviations 
outlined in Section 2. Overall, the system performed admirably, 
and the same dual -wedge system will be on board for the next 
flights beginning in the summer of 1984. The equations developed 
here for positioning and alignment will be used for these tests 
and will be evaluated in a ground-based experiment in early 1984. 
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